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We present the first evidence of charmless decays of the B° meson, the decay B°^</>0, and a mea- 
surement of the Branching Ratio BR(_Bg— *4'4>) using 180 pb" 1 of data collected by the CDF II exper- 
iment at the Fermilab Tevatron collider. In addition, the BR and direct CP asymmetry for the B + ^ 
<f)K + decay are measured. We obtain: BR(B°-><^>) = (14 t%(stat.) ± 6(syst.)) x 10~ 6 , BR(£+^ 
4>K + ) = (7.6±1.3{stat.)±0.6(syst.))xl0- 6 , and A C p(B + ^<t>K + ) = -0.07±0.17(stat) t ° 3 2 { s y st -)- 
Both decays are governed in the Standard Model by second order (penguin) boss's amplitudes. 

PACS numbers: 13.25.Hw 14.40.Nd 



Recent measurements |l|, |2| of CP asymmetries for 
B°^(j)Kg and B°^rj'Kg suggest deviations from Stan- 
dard Model (SM) predictions in B meson decays involv- 



ing b—>sss transitions, which also generate both B^—x/xp 
and B + — ► chK + decays. Charmless decays of B° and 
B + mesons provide among the most stringent tests of 
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the CKM matrix. Charmless B® decays, not yet experi- 
mentally established, offer important additional ways of 
testing the theory. Of particular interest for the current 
and next generation experiments at hadron colliders are 
charmless P" to Vector- Vector decays into self-conjugate 
final states, like P° — > cjxf). Exploiting angular correlations 
between the final state particles in B° s — > VV decays al- 
lows the statistical separation of the CP-even and CP- 
odd components of the decay amplitude. With sufficient 
statistics, this decompostion allows the measurement of 
the B® decay width difference (Ar s ), CKM studies, and 
tests of decay polarization predictions In particu- 
lar, the P° — > <jxj) decay, which proceeds through a pure 
b^sss, has been considered as a probe for New Physics 
scenarios (e.g. in b — ► s penguin dominated de- 
cays) recently advocated to explain the observed b-^sss 
asymmetries. The measurement of the branching ratio 
BR(P"— xjxf) gives insight into the size of penguin ampli- 
tudes as well as tightens the constraints on poorly known 
form factors. Recent calculations predict BR(P^ — > <fr<f>) 
between 10 x 1(T 6 5] and 37 x 1(T 6 0. Within the SM 
the B + ~^<pK + direct CP asymmetry (A CP ) is predicted 
not to exceed a few percent 0. No single experiment Q 
has yet reached the sensitivity to detect this asymmetry 
Further measurements of decay rates and Acp may help 
identifying the origin of the observed deviations from SM 
predictions in b— >sss penguin dominated B decays. 

We report the first evidence of P° — > (f>4> decays, and 
present the first measurements of the CP-averaged BR 
and Acp for B + ^<f>K + 9] at hadron colliders. To can- 
cel the uncertainty in the B hadron production cross sec- 
tion and to reduce systematic uncertainties on detector 
efficiencies, the branching fractions are extracted from 
ratios of the decay rates of interest normalized to the 
established B°^J/ip(/> and B + ^ J /^K + decay modes, 
which are characterized by the same number of decay 
vertices and charged tracks in the final state as in the 
signal modes. 

In this analysis we use 180 pb _1 of pp collision data 
at y/s — 1.96 TeV collected by the upgraded Collider 
Detector (CDF II) at the Fermilab Tevatron. The com- 
ponents of the CDF II detector pertinent to this analy- 
sis are briefly described. A more complete description 
can be found elsewhere lid. We use tracks in the pseu- 
dorapidity range \rj\ < 1 [ljj reconstructed by a silicon 
microstrip vertex detector (SVX II) and the Central 
Outer Tracker (COT) which are immersed in a 1.4 T 
solenoidal magnetic field. The SVX II detector consists 
of double-sided sensors arranged in five cylindrical lay- 
ers. Surrounding the SVX II is the COT, an open cell 
drift chamber with 96 sense wires. The integrated charge 
collected by each wire provides a measurement of the spe- 
cific ionization (dE/dx) for charged particles, allowing a 
separation equivalent to 1.4 Gaussian a between 7r and K 
for pt > 2 GeV/c. A set of planar drift chambers, located 
outside the calorimeters and additional steel absorbers, 



is used to detect muons within \rj\ < 1 with high purity. 

A sample enriched with heavy flavor particles is se- 
lected by the three-level displaced track trigger. At 
Level 1, charged tracks are reconstructed in the COT 
by the eXtremely Fast Tracker (XFT) 0- The trig- 
ger requires two oppositely charged tracks with trans- 
verse momenta pr > 2 GeV/c and the scalar sum pri + 
PT2 > 5.5 GeV/c. At Level 2, the Silicon Vertex Tracker 
(SVT) 0| associates SVX II r-<f> position measurements 
with XFT tracks, providing a precise measurement of 
the track impact parameter (do), the distance of clos- 
est approach of the track trajectory to the beam axis 
in the transverse plane. Decays of heavy flavor parti- 
cles are identified by requiring two tracks with 120 < 
do < 1.0 mm and an opening angle 2° < \A<f>\ < 90°. 
A requirement L xy > 200 /im is also applied, where the 
two-dimensional decay length, Pxy, is calculated as the 
transverse distance from the beam axis to the two track 
intersection projected onto the total transverse momen- 
tum of the track pair. A complete event reconstruction 
is performed at Level 3, where the Level 1 and Level 2 
trigger requirements are confirmed. 

B candidates are reconstructed by detecting <j> — > 
K + K~ and J/ip — > decays. In the latter case 

at least one of the muons has to be identified in the 
muon detectors to suppress contamination from other 
two-body J/tp decays. At least one pair of tracks ("trig- 
ger tracks"), has to satisfy the trigger requirements. B + 
(B®) candidates are formed by fitting three (four) tracks 
with pt > 0.4GeV/c to a common vertex. Requiring a 
good vertex fit x 2 reduces background from mismeasured 
tracks. Combinatoric background is reduced by exploit- 
ing several variables sensitive to the long lifetime and 
relatively hard px spectrum of B mesons and the isola- 
tion of B hadrons inside 6-quark jets. For this purpose 
we define the quantity Ip as the ratio of the B + candi- 
date pr over the total transverse momenta of all tracks 
within a cone of radius R = \J A77 2 + A<p 2 — 1 around 
the B flight direction. Requiring the B flight direction to 
extrapolate back to the beam axis decreases background 
from partially reconstructed decays. The cut values on 
the discriminating variables are optimized by maximiz- 
ing S I \/S + B for the already observed B + ~^> (f>K + signal 
and 5/(1.5 + V~B) for P^— >0</> whose branching ratio is 
unknown. The latter choice is equivalent to maximiz- 
ing the potential to reach a 3tx observation of a new sig- 
nal 0]. The signal (<S) is derived from a Monte Carlo 
(MC) simulation [T3| of the CDF II detector and trigger 
that uses the B meson momentum and rapidity distribu- 
tions from Ref. which were matched to CDF data. 
The background (B) is represented by appropriately nor- 
malized data selected with the same requirements as the 
signal except for the two-kaon invariant mass lying in the 
4> sideband region: 1.04 < tukk < 1-06 GeV/c 2 . 

We discuss first the B + — ► tfiK + analysis. The op- 
timization results in the following requirements: ver- 
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FIG. 1: rriKKK (left) and uikk (right) distributions for B + ~ -> 
4>K + with projections of the total likelihood function and its 
components. B + — > K + K~ K + is the sum of non-resonant 



and B + -> / (980)A: + contributions, while B" 1 
includes non-resonant and B + — >_ft'*°(892)7r + . 



tex x 2 < 8, B + decay length L xy > 350 /im, B + 
reconstructed impact parameter < 100 fim, isola- 
tion Iji > 0.5, non-trigger track transverse momentum 
> 1.3 GeV/c and impact parameter d^"^ > 120 /mi. 
The B^^>$K ± total yield and Aqp, defined as 



A 



CP 



N{B~ -> <pK~) - N(B^ 



N(B~ -► </>A-) + A(5+ -> 0A+) 



are extracted simultaneously from an extended unbinned 
maximum likelihood fit on four variables to the com- 
bined B + and B~ sample: the three-kaon invariant mass 
(ttikkk), the invariant mass of the <p candidate (mux), 
the cosine of the <j> meson helicity angle (defined as the an- 
gle between the K + momentum in the parent tf> rest frame 
and the momentum of the <fr in the B rest frame) and the 
measured dE/dx deviation from the expected value for 
pions for the lowest momentum trigger track. The like- 
lihood function has seven components: signal, partially 
reconstructed b— > <pX decays, combinatoric background, 
B+->K*°(892)ir+, with K*°-+K+ir-, B+^f (980)K+, 
with fo^K + K~, non-resonant B + ^K + K~K + and 
non-resonant B + — >K + tt~it + . The normalizations of the 
last three components are fixed to the B + ^K*°(892)ir + 
yield, determined in the fit, through their relative decay 
rates and detection efficiencies. For each component the 
likelihood function is the product of four one-dimensional 
probability density functions (PDF) of the fit variables, 
which are assumed to be uncorrelated. The ttikkk and 
nT>KK distributions are shown in Fig. ^ with projections 
for the different components. 

A combination of Monte Carlo simulation and side- 
band data is used to derive the PDF in each variable 
for the various fit components. For ttikkk the fully re- 
constructed signals are modeled by Gaussian functions. 
The mass and width of the B signals are determined 



from the fit in the case of KKK final states, while for 
A7T7T they are fixed to the values predicted by the simu- 
lation. We derive a parameterization from simulation for 
the partially reconstructed decays that populate the low 
mass side of the ttikkk distribution. An exponential 
plus a constant describe the combinatoric background. 
In the PDF for ttikk, the <\> resonance is described by 
a Breit-Wigner convoluted with a Gaussian resolution 
function, while the combinatoric background is modeled 
by an empirical phase space function. Shapes for other 
backgrounds are derived from simulation. The <p helic- 
ity PDFs for B decays are derived from simulation, while 
the combinatoric background PDF is modeled using data 
from the <j> sideband. The dE/dx PDFs for kaons and pi- 
ons are derived from a high statistics D°^> K~tt + sample 
obtained from D* ± decays. We find N^k = 47.0 ± 8.4, 
A C p = -0.07 ±0.17 and N K *o w + = 7.8 ±6.0 (statistical 
errors only) from which we estimate a B + — > /o(980)A + 
contamination of 11% under the signal peak. 

Candidates for the normalization mode, B + ~^ J/ij}K + , 
are selected with the same requirements as the B + — > 
(f)K + candidates except for the invariant mass of the two 
muons being within 100 MeV/c 2 of the J/tp mass ||. 
Using an extended likelihood fit of the m w jf and 
77i w distributions, we obtain a total yield of N^k — 
439 ± 22 (statistical error only). The asymmetry 
A C p{B+^J/ipK+) = 0.046 ± 0.050 is consistent with 
zero as expected for this mode. 

The relative B + ^c/>K + decay rate is calculated using: 



BR(B+^(j)K+) _ N^k BR {J/jj -> rifi) e^ K 



£<t>K 



e V>if ' 



BR (B+—>J/ipK+) N^ K BR (<j> -> K K) 

where e^if/e^K — 0.721 ± 0.011 is the ratio of the com- 
bined trigger and selection efficiencies derived from MC 
with a correction of about 5% due to the different trigger 
efficiency for muons and kaons as measured in unbiased 
samples. The efficiency for identifying at least one of the 
decay muons, e^ K — 0.81 ± 0.02, is obtained by weight- 
ing the expected px spectra in our signal with the single 
muon identification efficiency measured as a function of 
Pt in a sample of inclusive J/ip — > p + ii~ decays. The 
relative BR(B + ~^4>K + ) and Aq<p results are reported in 
Table |U Systematic uncertainties on signal yield (±1.4 
events) and asymmetry (±0.03, -0.02) are evaluated by 
varying the PDFs used in the likelihood fit, including a 
variation of the / (980) width from 40 to 100 MeV/c 2 . 
The uncertainties on the determination of efficiency and 
muon identification introduce a 5.6% relative error in 
the measurement of the branching ratio and have no ef- 
fect for Aqp- The uncertainties from BR((f>^rKK) and 
BR(J/Y— > fxfj,) contribute an additional 2% to the total 
systematic. We determine a ±0.005 uncertainty on the 
Acp measurement arising from a possible tracking effi- 
ciency asymmetry for K + and K~ as measured in mini- 
mum bias data. 

A "blind" search for B° B — ► <jxj> decays was performed 
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FIG. 2: Invariant mass distribution for Bg^cjxj) (left) and Bs^J/ip(f> (right) candidates. The arrows indicate the signal region 
for the Bj->^ search. The fit described in the text is superimposed on the m^KK distribution. 



TABLE I: Results for B + ^(j>K + and B°^cjxj> decays. World 
average 8] cj>— > K + K~ and J/ijj—* /i + /i~ BR have been used 
to obtain the BR ratio BR(B+ -» (j>K+)/BR{B+^ J/ipK+) 
and BR(B°->#)/BR(.B°-> J/ ipcj>). The first uncertainty is 
statistical, the second is systematic. 





B+^(j)K + 






Yield 


47.0 ± 8.4 ± 1.4 


7.3 ±2 5 


±0.4 


BR ratio 


(7.6 ± 1.3 ±0.6) x 10" 3 


(10 ±l±i 


) x 10~ 3 


A C p 


-0.07±0.17to;^ 







fixing the selection requirements and evaluating the com- 
binatoric background from independent samples before 
examining the signal region in the data. The signal is 
selected requiring two pairs of kaons having an invariant 
mass within 15 MeV/c 2 of the world average <f> mass 0. 
The optimized selection criteria are the following: ver- 
tex x 2 < 10, Bg decay length L xy > 350 /xm, £?° recon- 
structed impact parameter < 80 /xm, minimum <j> me- 
son transverse momentum pt, > 2.5 GeV/c and minimum 
of the two kaons' impact parameters from each 4> meson 
candidate (JqJ,,^ > 40 /im, d,Q 2 min > 110 /im, where 4>\ is 
the lower momentum <p candidate. The B® — > (fxfi candi- 
date mass distribution is shown in Fig. El In a region of 
±72MeV/c 2 around the world average |8j B° mass, cor- 
responding to a window three times the expected mass 
resolution, we observe 8 events. 

Two sources of background are expected in the B® sig- 
nal region: combinatoric background and B — > ipK* 
decays with the pion from the K*° decay treated as a 
kaon. The combinatoric background generates a smooth 
distribution in the invariant mass region close to mg«. 
Its contribution in the signal region is 0.35 ± 0.37 events 
(combining statistical and systematic uncertainties) , esti- 
mated using a background enriched sample where both 4> 
meson candidates have invariant mass lying in the cj> mass 



sideband region. The B ^ <fiK*° background results in 
an approximately Gaussian distribution underneath the 
B® signal. Its contribution, derived from simulation, is 
0.37±0.18 events where the error includes both statistical 
and systematic uncertainties, resulting in a signal yield 
of 7.3 5 events. The probability of a Poisson fluctua- 
tion of the background to the observed or higher number 
of events is 1.3 x 10~ 6 , corresponding to 4.7 a one-sided 
Gaussian significance. 

For the determination of BR(_B^ — > <j)(fi), a normaliza- 
tion sample of B® — > J/ip(j> decays is selected requiring 
one pair of kaons and one pair of muons within 15 and 
50 MeV/c 2 of the world average <f> and J/tp mass respec- 
tively. The other kinematic selection criteria are sim- 
ilar to the B® — > <f>4> mode. To extract the number of 
B® — > J/ip^ events, the candidate invariant mass dis- 
tribution, shown in Fig. El is fit with a binned maxi- 
mum likelihood function using a Gaussian for the signal 
and an exponential for the background. The fit returns 
Nm — 69 ± 10(stat.) ± b(syst.) events, where the sys- 
tematic error is evaluated using alternative background 
models for the low mass region where partially recon- 
structed B decays are expected. We subtract 3.7 ± 1.7 
background events from B° — > J/ijjK*° decays, with the 
pion treated as a kaon, as estimated from simulation. 

The B®—y(f)(fr decay rate is derived from the relation: 

BR (fl°->#) ^ BR ( J/V> -> wx) „ 
BR (B°->J/V><£) BR (cf> -> KK) e u ^ ' 

where e^^/e^^ = 0.821 ±0.015 is derived from simulation 
as in the B + ^ <f>K + case and , = 0.92 ± 0.05 is ob- 
tained using the pr spectra of the B®^ J/ip<p signal and 
the muon identification efficiency curve discussed above. 

The uncertainty on the B® — > J/ ifj/fi yield and back- 
ground evaluation contribute 8% to the relative system- 
atic error. The efficiencies significantly depend on both 
the polarization of the decay vector particles and the 
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assumed value of AT S since they affect the impact pa- 
rameter of the decay products on which the trigger op- 
erates. We vary the longitudinal polarization of the 
B° s -> (j)4> decay from 20% to 80%, the B a s -> J/ip<f> po- 
larization amplitudes within la of the CDF measure- 
ment and AT S in the range 0.06 < Ar s /r s < 
0.18 to assign a relative systematic uncertainty of 4%. 
Summing in quadrature all contributions, we estimate 
a total relative systematic uncertainty on the ratio of 
BR{B° -> <f>4>) to BR{B° S -► J/4>(j)) of 11%. To con- 
vert the ratio BR(S° ^)/BR(S° ■//■#) reported 
in Table |U in a measurement of BR(B° — * tjxf>), we use 
BR{B° S ^ J/ip<f)) = (1.38 ±0.49) x 10~ 3 , obtained from 
correcting the CDF measurement for the current 
value of f a /fd, the £>° to B° production ratio 0. We 
derive: BR(B° S -> 00) = (14 j^(siat) ± 6(syst.)) * 10~ 6 , 
where the systematic uncertainty includes a 36% contri- 
bution due to the uncertainty on BR(S°— > J/ij}(j)). 

In summary, we have used pp collision data collected 
with the displaced track trigger of the CDF II detector to 
study two fully reconstructed b^sss penguin dominated 
B decays. Using the world average 8] BR(B + ^J/^K + ) 
we measure: BR(B+ -> 4>K+) = (7.6 ± l.3{stat.) ± 
0.6(syst.)) x 10- 6 and A CP {B+ -> = -0.07 ± 

0.17(siai.) ln'o2( s ?/ s ^) which agree with previous mea- 
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surements |8JJ. We find the first evidence of a charmless 
Vector- Vector B° s decay and measure BR(£?° — > (fxj)) in 
agreement with the estimate of Ref. [j| and the recently 
amended calculation in |(J. 
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